Time-dependent unitary perturbation theory 
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We apply a time-dependent perturbation theory based on unitary transformations combined with 
averaging techniques, on molecular orientation dynamics by ultrashort pulses. We test the validity 
and the accuracy of this approach on LiCl described within a rigid-rotor model and find that it is 
more accurate than other approximations. Furthermore, it is shown that a noticeable orientation 
can be achieved for experimentally standard short laser pulses of zero time average. In this case, 
we determine the dynamically relevant parameters by using the perturbative propagator, that is 
derived from this scheme, and we investigate the temperature effects on the molecular orientation 
dynamics. 
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O ■ I. INTRODUCTION 

o 

Molecular alignment and orientation induced by an intense laser field continue to be a challenge to both experiment 
and theory. On the experimental side, it has already been shown that these processes have a large variety of applica- 
II 1 tions extending from chemical reactivity to nanoscale design 0,|2)0)SIEIE0I1|- From the theoretical point of view, 
several basic mechanisms aiming at limi ting the an gula r range of the molecular rotational motion have been derived 
in order to improve their control El llfl [ll [ Il2l Il3l Il4 | . In this context, one of the most efficient mechanism is the 
" kick mechanism" 0, UtI fl8L ITsi l20l |21| , which consists in a sudden angular momentum transfer to the molecule 
H ' by a so-called half-cycle pulse, which features a large asymmetry in the magnitude of the positive and negative peak 
values. The positive part, which is of interest in the process, is of short duration as compared to the free rotational 
• <-h , period of the molecule, and the negative par t is a weak and long tail that has a limited effect on the dynamics and 
can be neglected as a first approximation |17|. It is interesting to note that, no matter the intensity of such a field, a 
small parameter can be introduced by rescaling the time in the Schrodinger equation and that perturbative methods 
can therefore be applied with respect to this parameter |22j |. 

In other respects, it has recently been shown that time-dependent problems generated by short pulses can be treated 
by a time-dependent unitary perturbation theory (TDUPT ) Igp, 123. Igdl . which is the time-dependent version of the 
perturbation theory using averaging techniques [25l |2(J l28t l29l l3Ct l3l| . Analogous in its spirit to the interaction 
picture, this theory consists of a series of unitary transformations, which are aimed at rewriting the evolution operator 
as a product of other propagators. Moreover, order by order (Van Vleck) and superexponcntial (Kolmogorov-Arnold- 
Moser) algorithms can be derived. A detailed calculation scheme is provided in recent papers |2^,|2^| and its efficiency 
has been shown for a two-level Hamiltonian. 

The main purpose of the present article is to extend this study and to test the efficiency of the procedure as applied 
to the orientation dynamics of a diatomic molecule driven by an ultra-short electromagnetic field as a prototype. We 
demonstrate that the first order propagator can reproduce the time evolution in satisfactory agreement with exact 
results. We also show that an appropriate choice for a free parameter available in the TDUPT can help improving 
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the accuracy of the approximate propagator. Finally, we use the resulting perturbative propagator to investigate the 
postpulse molecular dynamics and to reveal its principal features. In this paper, we consider short pulses with no 
asymmetry in their temporal shape, which means that the time average of the electro mag netic field over this short 
duration is zero. Such pulses are actually those which are experimentally achievable |33.l33j. As the standard sudden 
impact approximation |18( is not a good starting point in this case, we construct a perturbative propagator using 
TDUPT. This enables us to determine the relevant parameters which control the orientation dynamics and to show 
that noticeable orientation can be obtained in this case. We also investigate the robustness with the temperature of 
such a mechanism and we remark that an efficient orientation may be reached for temperatures not exceeding 5 K for 
the LiCl molecule (which represents ten J states mainly initially populated). 

The remainder of this article is organized as follows : We outline the principles of TDUPT in Sec. |HJ paying 
special attention to the flexibility of the method. Starting from the expression of the Hamiltonian describing the 
molecule LiCl (within a rigid-rotor approximation) interacting with a linearly polarized laser pulse, we derive several 
propagators from TDUPT and compare them to the exact (numerical) one in Sec. IIIII The connections to other 
perturbation methods, like the Magnus expansion or the sudden impact approximation, are also presented. Section 
IIVI is devoted to a discussion of the parameters which control the orientation dynamics for experimentally available 
short pulses. For a temperature T = K, noticeable orientation lasting over 1 ps (or approximately one tenth of 
the rotational period) is obtained in the case of a strictly zero time averaged pulse. The results are also presented 
at T — 5 K and in spite of temperature effects which tend to decrease the orientation, it is shown that an efficient 
orientation can be obtained, but for shorter durations as compared to those observed at T = K. 



II. TIME-DEPENDENT UNITARY PERTURBATION THEORY 

This section recalls the principles of TDUPT as constructed in [22I I23I ] . The presentation followed here highlights 
the analogy of this procedure with the spirit of the interaction picture. We restrict ourselves to the first order, which 
we shall consider explicitly in the subsequent sections. To first order, it is worth noting that the general scheme 
described in this section is common to a time-dependent version of the Van Vleck and the time-dependent KAM 
perturbation procedures. The two methods differ however at higher orders. The reader is referred to refs. [22Ll23| for 
details. 



A. Description of the perturbation procedure 

Let H (t) be a time-dependent Hamiltonian, and Un(t, U) the corresponding evolution operator. Using atomic units 
(h = 1), Unit, U) is the solution of the Schrodinger equation 

d 

i^-U H {t,ti) = H{t)U H (t,ti), (1) 
at 

with the initial condition 

U H (ti,ti) = l, (2) 

where 1 is the identity operator. In order to solve Eq. by a perturbative scheme, we assume that a small parameter 
£ can be introduced in the following decomposition of the Hamiltonian H 

H(t)=H (t)+eV 1 (t), (3) 

where Hg is the unperturbed Hamiltonian and V\ the perturbation. We further assume that the evolution operator 
Uh of the time-dependent Hamiltonian Hq is known. 

The first step of the procedure consists in rewriting the propagator Unit, ii) in the interaction picture. Actually, one 
can introduce an additional parameter t p (the standard interaction picture being obtained for t p — 0) |!8ll34| 

U H (t,ti) = U Ho (t,t p )U Hl {t,ti;t p )U Ho (t p ,ti). (4) 

It turns out that the truncation of TDUPT at any order is strictly independent of t p [2i|. Hence, throughout the 
paper we set t p = U. In the interaction picture, Eq. Q reads 

i^-U Hl (t,ti) = eH^tpHd^U), (5) 
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where the Hamiltonian Hi is defined as 

Hi(t) =U\ Io {t,U)V 1 {t)U Ho {t,t i ). (6) 

Generally, the Hamiltonian H\ cannot in turn be partionned in a form similar to Eq. (j3J) with higher powers of e 
preventing thus the iteration of the procedure. However, a time-dependent version of unitary perturbation theory can 
be derived allowing the iteration through a series of unitary transformations. 

Following this scheme, the next step of the iteration consists in finding a unitary transformation : Ti(t;t 2 ) — 
e -ieH , i(t:t 2 ) ; where \y 1 } s self-adjoint, such that the resulting Hamiltonian H%, which will be explicitly defined below, 
can be decomposed in the form 

eH 1 (t)=sDi(t)+e 2 V 2 {t). (7) 

In this expression, Di is such that its evolution operator Ud 1 can be easily calculated, and e 2 V 2 contains no terms of 
order lower than 2 in e. Application of T\ to the evolution operator Uun according to the relation 

U Hl (t,U) = Ti(t;t 2 )U Ai (t,U;t 2 )Tl(U;t 2 ), (8) 

leads to the following expression for the Schrodinger equation [Eq. J5J] 

i^-Uh^U) = [ee iEW ^H x (t)e- isW ^ + i ^—— e - isW ^]U 6l (t,t i ). (9) 

For the sake of notation we omit the dependence on t 2 of T\, W\ and Ujj . The Hamiltonian Hi introduced in Eq. 
(JZJ) is accordingly defined as 

eHi(t) = ee isW ^Hi(t)e- isW ^ + i^— e ~ ieW ^ . (10) 

dt 

Expanding Eq. to first order in e, and taking into account Eq. (jJJ) one obtains the time-dependent generalization 
of the cohomological equation j2^| relating Wi to Di and Hi 

dWi 

—l = Hi{t)-Di{t). (11) 

The general solution to this equation reads 

Wx(t;t 2 )= [ du[Hi{u) - Di{u)l (12) 
Jt 2 

where t 2 is a free parameter. The important point to realize here is that this solution is not unique : The perturbation 
procedure offers the flexibility to control what part of Hi should be retained in W\ by a judicious choice of Di [23|. 
We consider explicitly several such choices in the next section. 

To summarize, with the help of the unitary transformation T\ we have shown that the evolution operator U fj 
satisfies the following Schrodinger equation 

i§^ Sl (t,U) = {eDi(t) + e 2 V 2 {t)} U &l (t,U), (13) 

formally similar to Eq J5J, but with a Hamiltonian Hi [given by Eq. Q] decomposed as H [Eq. (JSJ]. This allows for 
the iteration of the first step of the procedure. Writing U jj in the new interaction picture, we obtain 

U Ai (t,U) = U Dl (t,t q )U H2 {t,U;t q )U Dl {t q ,U), (14) 

where t q is arbitrary, and set to ti as, again, the final result is independent of its value. Uh 2 is the solution of the 
following equation similar to Eq. (JSJ 

d 

i-^U H2 (t,U) = e 2 H 2 {t)U H2 {t,U) 

= e 2 u]j x (t,t q )V 2 (t)U Dl {t,t q )UHA^ti)' (15) 
Up to this point, the procedure involves no approximation as it is just a sequence of unitary transformations. The 
first order approximation consists in replacing the evolution operator Uh 2 by the identity since it is associated with 
an Hamiltonian of order e 2 and is therefore equivalent to U~h 2 = 1 + 0(e 2 ). Combining Eqs. (gj, JSJ and l(T3)l. we 
finally obtain 

U H {t,U) = U„ (t,U)Ti(t)U Dl (t,U)Tl(U) + 0{e 2 ). (16) 

Equation (|16|) together with the various choices of D\ considered in the next section will be applied for studying 
molecular orientation through Sec. IIIII 
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B. Choices for improving the accuracy 



As it has been mentioned in the previous section, the principal point one can play with to favor the convergence 
of the perturbation procedure is the control of the terms of Hi to be kept in W x through the choice of D x . Three 
possibilities are considered : 

i) For a finite duration, examination of Eq. \12\ shows that no term must necessarily be kept in D\ to ensure a 
finite operator W\. The simplest perturbative propagator is therefore obtained by taking D\ — 0. Referring to 
Eq. I]16p. one is thus left with an evolution operator, U M , of the form 



tM 



U M (t, U) = U Ho (*, U) exp \te J duHi (u) j . (17) 

In this case, we remark that the first order propagator U M is completely equivalent to the usual first order 
propagator obtained with the Magnus formula for the Hamiltonian Hi [3^. The general relation between 
Magnus expansion and TDUPT has been derived in [23| . 

ii) Although not mandatory, secular terms, noted Si(t), can be taken into account in the definition of the operator 
D\. In the next section, we will see through numerical tests the role of such terms in the convergence of the 
perturbative propagator. Using an averaging method for time-dependent Hamiltonians 13, &i can be 
defined by the following expression 
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S i=J™ 7p / duH x {u). (18) 

T^+oo 1 J t _ T 

If H\ is an operator constant in time up to a time ti and with an arbitrary uniformly bounded dependence on 
time for t > ti, then it has been shown [2^ that Si = Hi(ti). The choice Di — Si allows one to determine the 
corresponding propagator as 

U Dl (t,U) = e - fe "i (*<)(*-**). (19) 

iii) Owing to the reduction of the error upon addition of secular terms in the definition of the operator D x , the 
question that naturally arises concerns the reduction of the error upon inclusion of other terms. Keeping in 
mind that the operator Di can be a solution of our problem as long as one is able to calculate the propagator 
Udi , we see that this second step is, however, far from being as obvious as the first one. Nevertheless, another 
simple and efficient solution consists in choosing the operator Di such that 

Di = Hx(ti), (20) 

where ti is a free parameter. Moreover, since Di does not depend on time, the calculation of Udi turns out to 
be very simple. We shall see that ti plays a significant role which can enhance the accuracy of the procedure 
by several orders of magnitude. This possibility, already present at the first order considered here, stems from 
the fact that to preserve unitarity we do not expand the exponentials, retaining thereby terms of higher orders. 

Having determined the operator Di , we can calculate the generator Wi from Eq. (|12|l and apply the general scheme 
of the procedure. However, we shall point out an additional adjustment possibility of this approach which is of crucial 
importance with the intention of improving the accuracy of the approximate propagator. This possibility consists in 
writing the Hamiltonian Hi defined by Eq. © in the following form 

eHi = sDi +e(Hi - Di). (21) 

The first step of the procedure can then be applied to the Hamiltonian eHi , where Hq is replaced by eDi and V x by 
Hi — D x . In this interaction picture, one can rewrite the propagator Uhi as 

U Hl (t,ti) = UD^top^t^topD^toA), (22) 

where to is a free parameter (set to ti from now as the final result is ^-independent) and Hi the resulting Hamiltonian, 
which plays the role of the Hamiltonian Hi in the general procedure described in Sec. Ill Al However, as Di has 
already been taken into account in the previous interaction picture, it is important to realize that the new operator 
Di calculated from the Hamiltonian Hi is generally taken as 0. The propagator obtained at first order in e with 
D\ = is hereafter referred as U 1 . 



How these adjustments are best performed for studying the orientation dynamics of diatomic molecules driven by 
a pulsed laser will be the subject of the next section. 
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III. APPLICATION TO THE ORIENTATION DYNAMICS OF DIATOMIC MOLECULES DRIVEN BY 

A PULSED LASER 

This section is devoted to the application of TDUPT for studying the orientation dynamics of polar diatomic 
molecules driven by an electromagnetic field. 



A. Description of the model 



We consider a molecule described in a rigid-rotor approximation interacting with a linearly polarized laser pulse. 
The model Hamiltonian is taken to be 

H = BJ 2 - noE(t) cos 9, (23) 

where J 2 is the angular momentum operator, B the rotational constant, fig the permanent dipole moment (for the 
sake of simplicity, the polarizability is neglected) and E{t) the electromagnetic field amplitude. 9 is the angle between 
the direction of the rotor axis and the polarization vector. The values /j,q = 7.129 D and B = 0.70652 cm -1 (value 
at the Li-Cl equilibrium distance |36j) are chosen so as to reproduce, at least qualitatively, the principal features of 
the polar diatomic molecule LiCl 11711 . We also recall that in spherical coordinates 9 (polar angle) and <f> (azimuthal 
angle), J 2 stands for the operator |3J| 

~ ~^9£)9 Sm 89 ~ (sm9) 2 deb 2 ' (24) 

Due to cylindrical symmetry, the projection m of the total angular momentum j on the field polarization axis is a 
classical constant of motion or a good quantum number. We shall consider, in particular a pulse shape E(t) of the 
form 

m = (E Q sin 2 (7r|) sin(27r/|) if < t < S (25) 
^ ' \ elsewhere, 

where Eq is the peak amplitude of the pulse, f/S its frequency and 5 its duration. Note that this function has been 
commonly used in the literature to describe the kick mechanism 0, 0, . The considerations developed below are 
not restricted to that pulse shape. In Sec. IIVI we shall be dealing with zero time-averaged pulses, which correspond 
to an integer value of the parameter / entering Eq. I|25|) . Experimentally, achievable values for Eq and 5 are taken to 
be Eq = 1.5 • 10 6 V ■ cm" 1 and 8 = 1 ps. This field duration is furthermore smaller by one order of magnitude than 
typical molecular rotational periods (1/B ~ 10 ps). 



B. Preliminary calculations and derivation of the perturbative evolution operator 

The time-dependent Schrodinger equation describing this system is 

i^U H {t,U) = {B J 2 - fi E(t) cos9} U H (t,ti). (26) 

Its evolution is governed by two characteristic times : the rotational period T rot = 1/ \J3 ( J 2 )] of the free molecule 
(where ( J 2 ) is the mean value of the operator J 2 ) and the pulse duration 6. Moreover, if \tp(t)}, the wave function 
of the field-free system at time t > 5, is one of the spherical harmonics \j, m) then this period can be written in the 
form 

(27) 



We remark that T rot decreases with the value of the quantum number j. 

In the sudden limit, i.e., for short-pulse duration with respect to the rotational period, a small dimcnsionlcss parameter 
e can be introduced as 



s = B5, 



(28) 
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which amounts to be e ~ 0.1 with the numerical values of B and 8. From a practical point of view, provided the 
pulse duration be sufficiently small, it is expected that this formulation allows one to study the dynamics of the 
system even in the presence of large non-perturbative pulse areas. Such a conclusion true for a two-level system 
has to be considered more carefully here in the sense that higher peak amplitudes induce higher rotational population 
leading to shorter periods, such that the choice of 8 actually depends on the pulse area. 

Rescaling the time in the form r = t/S, we obtain for the Schrodinger equation l|2tj|) 

i-^U H (T,n) = {-E r (T)cOs9 + Ej 2 }U H (T,Ti), (29) 

where E r (r) = [Iq8E{8t). We note Eq t = (j,q8Eq the peak amplitude of E r with a typical value of E$ r ~ 30 and we 
introduce the following dimensionless time = and Tf = 1. The pulse duration is then Tj < r < r/. Here, it is 
also important to realize that the small parameter e is independent of the choice of Eo r such that the perturbative 
method in consideration is consistent with very strong fields. 

The procedure described in Sec. Ill Al is applied to the Hamiltonian involved in Eq. 129(1 . In this particular case, 
we can first define the operators Hq and V\ of Eq. © as 



H (t) = -E r (r) cos 9 
Vi — J 



2 (30) 



hence 

U Ho (T,Ti) = e iA ^ cos0 , (31) 

where A(t) = duE r (u) is the pulse area in the interval [r,-, t]. Application of the first step of the procedure [Eq. 
©] leads to the following expression for Hi 

Hi(t) = e -iA{r) cosO j2 e iA{r) cob 6_ (33) 

Using the Campbell-Hausdorf formula [37| 

e c Be-° =B+±- l [C,B] + j l [C,[C,B]]+-.., (33) 

and the commutation relations 

[J 2 , cos 9] = 2(cr e +cos6»), (34) 
[ere, cos 9} = cos 2 9 - 1, (35) 

where erg = sin one can rewrite Eq. I|32|l exactly as [3l| 

Hi(t) = J 2 + 2i A(T){a g + cos 9)+ A 2 (t){1- cos 2 9). (36) 

Finally, to first order in e, the propagator of Eq. I|l(j|) is written in the following form for Tj < r < Tf 

U H (r,n) = e *^( T ) cos ®T\{t)U £ d 1 (t, Tj)T^(rj) + 0(e 2 ), (37) 

The choice for the operators D\ and T\ will be explicited in the next section. For r > rt, the propagation is assumed 
to be free 

U H (r,r f )=e- is i—f) J \ (38) 
C. Numerical tests 

We are now in a position to check the accuracy of the perturbation scheme through the comparison of wave 
functions obtaine d by the various TDUPT propagators with those obtained from the accurate full numerical split- 
operator method |38l l39| . The choices described above for Di and for the decomposition of the Hamiltonian Hi , as 
well as the values of the time parameters T\ and r 2 , lead to different evolution operators, the merits of which are 
numerically checked hereafter. More precisely, for Tj < r < r/, we quote : 
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i) The choice D\ = which leads to the first order Magnus propagator through the use of Eqs. (|17|l . (|31|) and 



(39) 



ii) The choice D\ = Si with a secular term S\ [Eq. (|18|l ] that can easily be evaluated through Eq. l|3t)|) with the 
result in the case where A vanishes for r < Tj 

Si = J 2 , (40) 
leads to the propagator II s calculated using Eq. (|19|l 

U s (r,Ti) = e a W cose exp|-iey T du[JIi(u)- J 2 ] | 

xe-^tr-r,),/ 2 exp |_ . £ /" 2 du r Hi ^ _ j2j I ( 41 ) 



(42) 



hi) The decomposition defined by Eq. I|21|) with Z?i = i?i(n) as in Eq. H2U|) leading to 

x exp | -*ey due^"-^ 1 ^ [ff^u) - JZi( n )] e -«(«-«)^i(ri) J , 

In the particular case t\ = Ti, one deduces from Eq. (|36|) that 

H\{t\ — Ti) — J 2 . (43) 

It is worth mentioning that all these propagators are close to the Magnus one in a sense specified in |2^| . However, from 
a pratical point of view, the fact that the operators J 2 , ag and cos# do not commute leads to different propagators. 
We also recall the form of the propagator U SI that has already been used for this problem in the sudden impact 
approximation [3] • Rewriting the initial Hamiltonian [Eq. (|29|l ] in the interaction picture involving the operator J 2 
and neglecting the molecular rotational motion during the pulse, yields at time r/ 

U SI (T f , Tl ) = e -Hr f -T h )J 2 e iA(Tf) C o S d e -ie(T h -ri)J\ (44) 

where is a parameter. Numerical tests show that the best choice for this free parameter is = (tj — t.;)/2. It 
is interesting to note that this parameter t^ although introduced through an interaction representation (see Eq. 
(@J and remark below this equation) does affect the final result. This is due to the fact that in 18] an additional 
approximation is made which assumes an impulsive character for the perturbation. 



In order to measure the accuracy of these propagators, we define the error A at the end of the pulse as 

A=||^r / )-^(r / )|| 2 (45) 

where ip ex is the exact wave function computed by solving the time-dependent Schrodingcr equation with the split- 
operator method and ^ the one obtained by applying the propagator under consideration. In the numerical cases 
studied below, the initial wave function V'('ri) is taken as the ground rotational state \j = 0, m = 0) of the molecule. 
Note that this state can be experimentally prepared, for instance, by laser cooling methods |40| . 

The remarkable accuracy of TDUPT propagators is clearly illustrated in Fig. n an d|21 which display the logarithm 
of A as a function of the parameter e for the following evolution operators U M , U SI , II s and U 1 . Numerical values 
are Eq t — 1 and / = 0.5 in Fig. ^and Eo r = 20 and / = 2 in Fig. [21 The pulse area is zero when / is an integer, 
which allows us to consider stronger fields in this case. It is also noted that, for an half-cycle pulse (/ = 0.5), the 
maximum peak amplitude that can be experimentally achieved is smaller than Eq = 1.5 • 10 5 V ■ cm -1 which 
corresponds to a parameter E$ r of the order of 3. As could be expected, the general trend is that the error decreases 
as the parameter e decreases. In this way, the exact result is reproduced within an accurac y b etter than 10 -3 when 
e < 0.5. More unexpectedly, the agreement is still quite good for larger values of e (see [23| for an explanation). 
Indeed, the error for e = 1 is as small as A ~ 10~ 4 for U 1 and / = 2, which is particularly impressive considering that 
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only one iteration is used. Moreover, we remark that the secular terms do not decrease the error of the perturbative 
propagator II s in comparison with U M whereas, with the addition of the last adjustement [Eq. l|21[l ]. A is roughly 
reduced by at least, two orders of magnitude for U 1 . For e < 1, U 1 is much more accurate than the sudden impact 
evolution operator U SI . This difference is particularly striking when the pulse area is zero (/ = 2, Fig. |2J). We next 
analyze the logarithm of A as a function of the parameter E 0r for the propagator U 1 . Such a plot appears in Fig. |3 
where numerical values for e and / are taken to be e = 1 and f — 2. As could be expected, the most salient feature 
of Fig. [21 is the fact that the error increases with E Qr . This means that the rotational population is strongly modified 
and that higher J's have to be taken into account. 

The second part of this section concerns the role of the free time parameters t\ and ti in the lowering of the error. 
First notice that U M and U 1 do not depend on t^. For T2-dependent TDUPT propagators as U s , it can be shown 
that this dependence is not dramatic and the best choice is T2 = r-j [2^ |. 

Figure Q] displays the logarithm of A for the propagator U 1 as a function of the time parameter n . Numerical 
values for e, / and E$ r are taken to be e = 1, / = 2 and Eq t = 10. In this example, it is clear that the best choice 
for this parameter is T\ ~ 0.11 (or t± ~ 0.89). This degree of freedom can be optimized without any, a priori, knowl- 
edge of the exact solution by locating the minimum with respect to t\ of the eigenvalues of an appropriate operator |'2.'"{| . 
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FIG. 2: Same as Fig[T] but for / = 2 and E 0r = 20. 



As explained in this section, several choices and attempts can be made to improve the TDUPT propagators. 
Nevertheless, the best choice of the operator D\ or the best set of time parameters remains dependent upon the 
context. As long as no dynamical information is lost, the best choice is the one for which the perturbative propagator 
is as simple as possible. 



IV. ORIENTATION DYNAMICS BY SHORT PULSE OF ZERO TIME-AVERAGE 



A. Zero rotational temperature 



We first consider the limit of low rotational temperatures. Having demonstrated the validity and the accuracy of 
the perturbative propagator, we now use it to investigate molecular orientation dynamics after the pulse is over. More 
precisely, we are looking for the relevant parameters of the laser imp roving the orientation of the molecule. The mean 
value (cos#) r is taken as a quantitative measure of orientation |171 141| 

(cos0) T = (<Kt)|cos0|V>(t)) • (46) 

It is to be noted that a good orientation is obtained for large absolute values of (cos8) T and that the measure used 
does not take into account the temperature effects jljj, |2(J, |42| , which will be investigated in the second part of this 
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FIG. 3: Common logarithm of the error A as a function of the adimensional 
parameter Eo r for the propagator U 1 in the case e = 1 and / = 2. The free time parameter ri is fixed to 0. 

section. 

In order to highlight the role of TDUPT in the understanding of molecular orientation, we consider short pulses 
with symmetrical temporal shape, which means that the time avera ge o f the radiative field over this short duration 
is zero as it has to be for a freely propagating electromagnetic pulse |32| . To our knowledge, it has never been shown 
that such pulses can be used to obtain a good orientation. Moreover, according to the sudden impact approximation, 
so far used in the literature to describe the kick mechanism [Tj, Ha. l2ll | , such (zero time averaged) pulses would lead 
to post-pulse dynamics without any orientation effects. Indeed, the sudden impact propagator U SI is given by the 
following expression [Eq. i|44p] for r > 17 

JjSIr T ^ T .\ = e -ie{T-T h )J 2 e %A{T f ) cosd e -ie{T h -Ti) J 2 ^ r^j\ 

If A(rf) — and ip{ T i) = \j> m )t it is clear that (cos6>) r = when the pulse is turned off. This is also seen in Fig. [S] 
which displays (cos#) r as a function of time r. 

Using the propagators we have just derived with TDUPT, we calculate (cos#) T after the pulse. We consider the 
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1 

FIG. 4: Common logarithm of the error A as a function of the adimensional 
time parameter ri for the propagator U 1 in the case e = 1, / = 2 and Eo r — 10. 

most accurate evolution operator U 1 for r > tj and take n = Tj, so that Eqs. H42[) and (|43[) yield 

^ (T , T() . exp {_ i£ £ [Hl{u) _ J2] .-^J . 

From Eq. one deduces that 

ffi(u) - J 2 = 2iA(r)(ae + cos6>) + A 2 (r)(l - cos 2 6>), 

where we recall that ag = sin 6d/d9. 

Let us define the following operator which appears in the second exponential of Eq. (|48|l 

F = i j*' due l£ ( u ~ Ti)j2 {H^u) - J 2 } e -^-r,)J\ 
We rewrite accordingly the propagator as 

U I (r,n) = e-^-^ j2 e- sF . 



(48) 



(49) 



(50) 



(51) 
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By virtue of the action of the operators J 2 , cos 9 and og in the basis of the spherical harmonics 



J 2 \h m ) = j(J + 1) \j,™>) , 
cos0|j,m) = c j+lm \j + l,m) +c jm \j ~l,m) , 

°e | j, m) = j'Cj+i m |i + 1, m) - (j + l)cj m \j - 1, TO) , 

1/2 

where Cj m = [(j — m)(j 4- m)/(2j — l)(2j + 1)] , it follows that F |j, to) is of the form 

F \j, m) = «j m |j - 1, to) + bj m \ j + 1, to) + aj TO | j - 2, to) + /3j m | j + 2, to) + 7j m | j, to) 
One readily obtains the coefficients 



bj rn 

Pj m 
Jjm 



T f 



= 2jc jm / duA{u)e- 2 ^ u '^\ 



-2(j + l)c j+lm J duA(u)e 2 ^ +1 ^ u -^\ 

"Tf 



ICj — lmCjin 
*(1 — c jm ~~ c j+lm 



duA 2 (u)e 



-2ie(2j-l)(u-Ti) 



Tf 

duA 2 (u). 



(52) 
(53) 



(54) 



For a system initially prepared in the state \j, m) 



\^{n)) = \j,m) , 



(55) 



we describe the orientation dynamics at time r with the mean value (cos#) T computed with the propagator U 1 
according to Eqs. l|4^|) and (|5T)l 



(costf) T = ^j,m 
Using the Campbell-Hausdorff formula, Eq. (|33|l , we obtain 



J,m). 



(56) 



(cosw; r = eij,m 



[F, 



cosf 



-is(T-T,)J 2 



j,m) + 0{e ) 



Aejc 2 



= Ae{j + l)c 2 m / duA(u) cos (2s[j + l][u - r + tj]) 



dwA(u) cos (2e?'[tt - r + r 4 ]) + 0(e 2 ). 



(57) 



Note that the higher order terms can be neglected as the propagator U 1 is constructed using the first order TDUPT. 
In other words, at this level of approximation we have already discarded terms of order e 2 . 

For very low rotational temperatures, the only rotational level initially populated being \j — 0,to = 0) 40], one 
deduces that 



(cosf 



e I cos(2e[r - n]) I duA{u) cos(2eu) + sin(2e[r - rj) / duA(u) sin(2eu) \ + 0{e 2 ). (58) 



Introducing the Fourier transform A of the pulse area 

1 f Tf 



A(k) 



'2tt J Ti 



duA(u)e 



we rewrite Eq. (|57|) in the form 



4V2^ 



(cos( 



A(2e) cos [2s[t - n] + argL4(2e)]) + 0(e 2 ). 



(59) 



(60) 
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FIG. 5: Mean value (cos6 l } T as a function of the adimensional time r computed with the split-operator method (solid line), 
the sudden impact propagator U SI (dotted line), and the perturbative propagator U 1 according to Eq. 1561 1 [dash-dotted line] 
and Eq. [dashed line]. Here T = OK for e = 0.5, / = 2, E 0r = 50 and the pulse is on for < r < 1. 

In Figure|3]we display the mean values (cos0) T computed with the propagator U 1 according to Eq. I|56|l and according 
to its first order contribution as given in Eq. H6()(l . The mean values obtained with the sudden impact propagator 
of Eq. (|47|l and purely numerically with the split-operator method are also depicted. It is seen that the first order 
expression of Eq. (|60|l brings a significant improvement with respect to the standard sudden impact approach which 
predicts no post-pulse orientation effect in the case of short pulses of zero time-average considered here. 

It is worth noting that 2e is the rescaled frequency between the first two rotational levels and that the leading 
revival structures [ija] are well described by this approximation up to large values of e, at least from a qualitative 
point of view. The other fundamental frequencies are neglected to first order in e, which physically means that the 
corresponding levels with higher j's are not strongly populated in comparison with the first two rotational states. 
From Eq. t|f>()p one deduces that the range of the post-pulse orientation mainly depends on e and |A(2e)| which are 
the two dynamically relevant parameters one can play with to control (cos#) r . 

Moreover, thanks to a careful choice of the laser parameters which will be discussed below, it is shown in Fig. 
[3] that a noticeable orientation can be achieved with zero time averaged pulses. Indeed, in this example the value 
| (cos 0) | ~ 0.5 is reached and lasts for a time larger than 1 ps. To our knowledge, it is the first time that such an 
orientation is observed in this case. 

We next analyze the choice of the different parameters. For e sufficiently small, |A(2e)| is well approximated by 
J Tf duA(u). It is instructive to consider the case where the pulse shape is given by Eq. (|25|l for which we see how 
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this latter quantity depends on the frequency and the peak amplitude 

rn [ E orlh( T f ~ T *) if/ = l 

/ duA(u) = I (61) 
Jr< ( E orh7U^i T f - T *) if / > L 

It follows from Eqs. (|60[1 and Ij61|l that the orientation decreases as the frequency / increases. In the high-frequency 
regime, where / ^> 1, we see that no orientation can be obtained as has already been shown using a high-frequency 
Floquet approach |4ij. On the other hand, in the range of validity of Eq. ijSUjl . it is important to realize that the 
increase of e and Eo r involves a better orientation which, however, lasts for shorter durations (the period of the 
motion is equal to tt/e). Experimentally, this point is of crucial importance because the possibility to perform, for 
instance, stereodynamically sensitive chemical reactions by using this orientation depends on this period of time. For 
practical purposes, this duration has to be larger than 1 or 2 ps. Finally, numerical values for e, E$ r and /, which 
fulfill these conditions, are taken to be e = 0.5, E 0r = 50 and / = 2. We stress that the conclusions on the choice of 
parameters are general, and not restricted to the particular case of molecule (LiCl) considered as an illustration. 



B. Finite rotational temperature 



The results presented above were based on the assumption that the rotational temperature T was zero. We now 
investigate the temperature effects on molecular orientation. In this case, a thermal average over the rotational levels 
has to be taken into account. The quantitative measure of the orientation is then given by [T^. I2H 142^ 

« cosfl > r = 1 gexp ( ~ B ^ + 1] ) E <cosfl) T , (62) 

3 m=—j 

where k B is the Boltzmann constant and Q the partition function 

Q = E(2i + l)expf -^_ +1] ). (63) 



The result for LiCl molecule under the effect of the previous pulse is shown in Fig. which displays the thermally 
averaged mean value <C cos 9 3> T as a function of the time r for the rotational temperature T — 5 K. 

We first note a decrease of the orientation with increasing temperature, this point has already been mentionncd 
in previous studies [T^. I20I l42| . Moreover, as for the study at T = OK, an approximate analytical formula can be 
derived for <C cos 9 ^> T . Using Eq. J57J|, straightforward calculations lead, to first order in e and for r > r/, to : 

8e + °° / Bi \ ( Bj 2 \ f Tf 

<cosfl> r =— ^jc. sinh^^Jexp^-^J J duA(u) cos (2ej[u - r + n\) + 0(e 2 ), (64) 

where c, is defined by 

+j 
m=—j 

Introducing as previously the Fourier transform A [Eq. Q59[l]. -C cos 9 3> T can then be rewritten in the form 

8V2^e . . u / Bj \ ( Bf 
< cos6> > T = — — 2^ JCj smh — — I exp 



Q ^ J J \k B T r V k B T 
A(2sj)\ cos (2ej[T - n] + arg[A(2ej)]) + 0{e 2 ). (66) 



The corresponding orientation dynamics, displayed in Fig. fairly reproduces the position and the value of the main 
extrema of <C cos 9 ^> T which is the objective of this analysis. Analysing along the same lines as T — K the role 
of laser parameters, one deduces that a noticeable orientation can be obtained for the following choice of numerical 



15 




FIG. 6: Thermally averaged mean value <C cos 9 2> T as a function of the adimensional time r for the rotational temperature 
T = 5K for e = 0.5, / = 2 and Eor ~ 70. The pulse is on for < r < 1. The solid line corresponds to the split-operator 
method and the dash-dotted line represents the result obtained from the perturbative propagator U 1 according to Eq. 1661 . 

values : e = 0.5, / = 2 and E$ r = 70. In this way, an orientation efficiency of <C cos# 3> r ~ 0.5 with a duration larger 
than At = 0.3 ps is achieved, which also shows the robustness with respect to temperature of this mechanism. At is 
the duration over which <C cos9 ^> T remains larger than 0.3. To our knowledge, this result corresponds to one of the 
best reported in the literature. On the other hand, apart from these extrema, the thermally averaged <C cos 6 ^> T is 
close to the value zero, which corresponds to no orientation. More precisely, as can be clearly seen in Eq. (|66|l . the loss 
of orientation between each maximum is due to the different periods n/sj (j ' > 1) of cos (9 3> T , corresponding to 
the different rotational frequencies of the molecule. Such an effect is basically expected since, for a freely propagating 
pulse, the time average of (cos#) r over a rotational period is zero 

1 f t+T tat 

— du{co S e) T (u)=Q. (67) 

^rot Jt 

V. SUMMARY 

This article has focused on the application of the TDUPT for studying orientation dynamics of diatomic molecules 
driven by pulsed laser fields. Numerical tests have demonstrated the efficiency of the proposed procedure. A basic 
feature and advantage of the method under consideration, where the small perturbation parameter is the short pulse 
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duration, is that it allows, through an analytical description of the molecular dynamics, for a thorough interpretation 
of intense laser induced orientation. 

Until now molecular orientation with short laser pulses had been only envisaged with non-zero time-average of the 
electric field as the standard sudden impact propagator predicts no post-pulse orientation when this time-average 
vanishes. Since a free propagating electromagnetic wave must posses a zero time-averaged electric field, half-cycle 
pulses have been used. They are composed of a short intense pulse with non-zero time-average (which induces the 
molecular orientation) , and a long weak tail (which is neglected) . In this paper we considered symmetrical short laser 
pulses with zero time-average which are easier to produce experimentally. We constructed a perturbative propagator 
that enables us to elucidate the post-pulse orientation dynamics. For zero time averaged short pulses and low rotational 
temperatures, we have shown that the orientation may be significant depending on two leading parameters. On the 
one hand there is the adimensional parameter e = B5 where B is the rotational constant and <5 the pulse duration, 
and on the other hand, A(2e), the Fourier transform of the pulse area evaluated at the rescaled frequency between 
the first two rotational levels. The orientation is proportional to |^4(2e)| which for a prototype of symmetrical pulses 
is itself proportional to the peak amplitude and inversely proportional to the frequency. Finally, we investigated the 
effect of temperature and showed that a good orientation with a shorter duration can be achieved for finite but low 
temperatures (e.g., up to 5 K in the case of LiCl) by an adequate choice of these pulse parameters. 
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